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We investigate the long distance contribution via charmed hadron loops in the processes r/c, XcO 
and he decaying to baryon anti-baryon pairs, which are supposed to be highly suppressed by the 
helicity selection rule as a consequence of the perturbative QCD framework. With an effective 
Lagrangian method, our estimation result indicates that such hadron loops play an important role 
in these hadronic decays. It is a further test of the evading mechanism for the helicity selection rule 
in charmonium baryon-antibaryon decays. 
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I. INTRODUCTION 



, The complexity of quantum chromo dynamics (QCD) remains unsolved in many questions in the charmonium mass 



- '—I 

X 



region. This is the regime that numerous contradictory results between perturbative calculations and experimental 



^ observations were found. In particular, the pQCD expected helicity selection rule has been found violated in many 
^ ^ ' exclusive decay processes. According to this selection rule, some charmonium decay channels are supposed to be highly 
suppressed, such as J/ip — > VP, rjc VV and Xci — > VV ^ where V and P denote vector and pseudoscalar meson 
\ respectively. However, all these decays are found to be rather important in experiments In Ref. [5|, Feldmann and 
Kroll argued that charmonium mesonic decays could be classified into two catalogues. One is controlled by pQCD with 
I— I conserved hadronic helicity while the other is characterized by the helicity-selection-rule violation. It was proposed 
i-C ■ that a soft mechanism through a light-quark Fock component would account for the violation of the hadronic helicity 
conservation in exclusive decays. Other theoretical solutions for the underlying dynamics have also been explored 
^ in the literature [5. ^,10:1 . and detailed review of some of the related questions in charmonium exclusive decays can be 
D ; found in Refs. [iH1[11Hi1. 

. In this work we are going to investigate another three processes, i.e., rjc,XcO,hc — >■ YY, where YY represent the 
= 1/2+ octet baryon-antibaryon pairs, i.e. pp, AA, EE, and SS. These channels are also supposed to be highly 



: suppressed according to the helicity selection rule [1J|. However, it seems that the available experimental data do 



^ ' not support such expectations at all [J]. Some attempts have been made in order to understand this contradiction 
i [l5l - l20l |. For instance, in the processes of rjc and Xcji-J = 0, 1, 2) decaying into pp, a quark-diquark model for proton is 
introduced as a mechanism for evading the helicity selection rule. However, the obtained branching ratios for rjc and 
' XcJ can not simultaneously agree with the experimental da ta fisl [l6| . Constituent quark mass corrections were also 
introduced to account for the processes rjc, XcO: he — > pp [itI llSj. where it turns out that the obtained branching ratio 



\l for rjc pp is still much smaller than the experimental data. There are also some other proposed mechanisms for 
, understanding these processes, such as the mixing between charmonium state and glueball and the quark pair 
^-"^ ' creation model [2^. Basically, it is believed that non-perturbative mechanisms should be relevant in this scenario. 

In Refs. [1, Q, it was proposed that intermediate meson loop (IML) transitions can serve as a soft mechanism in 
charmonium decays. Such a long-distance interaction can evade the Okubo-Zweig-Iizuka (OZI) rule and result in 
violation of the pQCD helicity selection rule. In the framework of effective Lagrangians for hadron interactions, we 
can then quantitatively study the evasion of the helicity selection rule in various processes. As a further check of 
d ' this mechanism, we are going to investigate the role of the charmed hadron loops in r]c,XcO,hc — >■ YY in this work. 
Additional evidence for such an underlying mechanism should allow us to gain more insights into the QCD strong 
interaction properties in the intermediate energy region. 

The rest of this article is arranged as follows: In Sec. II, we will depict the effective Lagrangian method, and some 
relevant formulas will be given. In Sec. HI, we will present the numerical results and discussions. The conclusions 
will be summarized in Sec. IV. 



* xhliu@ihep.ac.cn 
t zhaoq@ihep.ac.cn 



2 



II. LONG-DISTANCE CONTRIBUTION VIA CHARMED HADRON LOOPS 



In Figs. [TJ [21 E] and HI the intermediate charmed hadron loops which serve as a long-distance soft mechanism are 
illustrated by triangle diagrams. Since one charmed baryon will be present in the loop, these diagrams are somewhat 



9|. We will consider the exchanges of the ground state 
2l| . There are several points that should be clarified 



different from the intermediate meson loops studied in Refs. 

= 1/2+ charmed baryons that belong to SU(4) multiplets 
for Figs. [man and H 

1) Jcc represents the decaying cc meson, i.e. rjc, XcO, or he- 

2) D(^) and Dj-g) only denote the flavor contents but do not contain the spin quantum numbers, and they will depend 
on the decaying mesons. 

3) Different intermediate charmed hadrons will appear in the loops for ?7c,XcOi^c ^ YY. For the r]c decay, there 
are three situations for D(s) and D(s): (a) D(^s'jD*^y (b) D*^^D(^s), (c) D'^^s^D*^^y For the XcO decay there are two: (a) 
-D(s)-D(s), (b) D*^^D'^^y For the he decay there are three: (a) D(s)DL)^ (b) DJ^s^Dfs^, (c) D'^^^D*^y This is due to the 
adopted effective Lagrangians based on heavy quark symmetry [2 21 123| . 

We list the relevant effective Lagrangians as follows: 



Ci ^ igiTr[Pl^eH2ilt^Hii] + /i.e., 
C2 = ig2Tr[RccH2iY' diiHii] + h.c, 
where the spin multiplets for these four P-wave and two 5-wave charmonium states are expressed as 



^ cc 
Rcr. 



l+j! 

2 



1 

V2^ 



1 



1 - ^ 



(?A^7m " '7c 75) 



1 - f 



2 y 2 

The charmed and anti-charmed meson triplet read 



(1) 
(2) 

(3) 
(4) 

(5) 
(6) 



where V and V* denote the pseudoscalar and vector charmed meson fields respectively, i.e. I'*^*-' = 
(£,o(*) r,+(*),i^+(*)). For the meson-baryon interaction Lagrangians, we follow the forms that were adopted in 
Refs. 



^YaVY 
^YaVY 



igY^VYYc'-fsYV + h.c, 

I^YaVY 



9y^v^yYc i^'V 



2mAr 



a^'''^JD* Y + h.c 



(7) 
(8) 



where Yc, 'D^*\ and Y denote the charmed baryon, charmed meson, and the corresponding nucleon or hyperon, 
respectively. The relevant coupling constants will be discussed later. 

With the above effective Lagrangians, we can now calculate the transition amplitudes illustrated in Figs. [1] [21 El 
and[4l For these diagrams, we take the convention of the momenta as Jcc{p) — ^ '^{qi)^[(l2)\Yc{q)\ — >■ Y{pi)Y{p2)- 




(1) (2) (3) 

FIG. 1: Charmed hadron loop diagrams that describe the long-distance transitions in Jcc — > pp. 
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FIG. 2: Charmed hadron loop diagrams that describe the long-distance transitions in Jcc — >■ AA. 

^ , yi « — Tni+ yr—* E"*" 




(1) (2) (3) 

FIG. 3: Charmed hadron loop diagrams that describe the long-distance transitions in Jcc — >■ E'E"* 



A. r)c YY 



We first consider rjc — YY , where three different amplitudes corresponding to those three intermediate states will 
contribute, i.e. M.a, M.b, and M.c- These amplitudes can be written down explicitly as follows: 



X u{p2) ( 7m + ^-^^3~^'^f^''l2 ] (4 + myj'fsvipi) 



2m N 



9 2 ? 9 9 9 

— my qf — mj, q2 — mjy 



f d^q ( q^q^ 



1 



1 



1 



9 9 9 9 9 9 

q~ - mf^ q{ - m^, q^ - 



X u{p2) ( 7r + » II!^ ^ f^r^qi ) + "lyj ( 7a + i ^f^ ^ tJAag?" ) "wbl) 



2mN 



2mN 



(9) 



(10) 




FIG. 4: Charmed hadron loop diagrams that describe the long-distance transitions in Jcc ^ 
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X — 2 2 2 2 —^W)- (11) 

Since the masses of ijcXcO, and /ic are below the threshold of D^*^D'^*\ intermediate mesons and can 
not be on-shell simultaneously. We thus phenomenologically introduce a form factor T{q^) as has been done in Refs. 
0, [12: Hlj to take into account the off-shell effects, 

where qi = q, qi, q2- The cut-off energy is chosen as Ai = rrii + aKqco, Aqcd — 0.22 GeV and is the mass of the 
corresponding exchanged particle. The form factor is also necessary for killing the divergence of the loop integrals, 
although it will give rise to mo del- dependent aspects of the calculations. Qualitatively, this type of form factor will 
converge the integrals faster than a dipole form factor and can be well combined with a vertex coupling constant as 
a phenomenological account of the non-local coupling form factor. Some further discussions will be given later in the 
numerical analysis. 

In this work, we do not include contributions from exchanging other charmed hadrons with higher spin or orbital 
excitations. These resonances have relatively larger masses and their couplings so far are unknown. It should be 
reasonable to only consider the lowest partial wave states based on the argument of locally broken down quark-hadron 
duality 

B. xco YY 

Two charmed hadron loops will contribute to XcO YY , and the corresponding amplitudes are 

^ ' ^ H<1% (13) 



9 9 9 9 9 9 

q^ - mf^ ql - q^ - 



... 2 / d q ( ^£g£\ / ^„ ql^q^ 



(2^)4 •'^M- V ^IJ \ rnl. 



X U{p2) [ 7a + i ga/3g2 j H + mY,) p + i gprglj v{pi) 

^ ^ Hq^), (14) 



9 2 2 2 2 2 

q-^ — my q^ — m^, q^ — mij~, 
where the form factor has the same form as that in rjc — > YY . 



C. /le YY 



Similarly, the amplitudes for he — >■ YY from those three contributing loops can be written down as follows, 



Ma = gh^vvgY^VYgr^VYertip) / TTTTT -.9^'' + ^2^ 

I^Y^VY 



1 1 1 2n 

2 _-Z2- „2 „ ^2_ „2 _ ^2_ -^(9 )' (15) 



Mb = gh^vvgY^vYgY^vY£-n{p) J 



4 ^ ' ml. 



X w(p2)75(^ + ™yj ( 7m + » ^2m^^ c^pt^gi ) '^(Pi) 
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9 9 9 9 9 9 



^ J-(g'), (16) 



ml,. 



X „ ^ 2 2 ^ 2 -^(g^)- (17) 

where e,j(p) is the polarization vector for and again the form factor has the same form as that in 77c YY . 

III. NUMERICAL RESULTS AND DISCUSSION 

Proceeding to the numerical results, we first discuss the determination of coupling constants. For the couplings of 
the charmonium states to charmed mesons, the expansion of the effective Lagrangians Ci and £2 gives the following 
relations in the heavy quark limit: 

m-D* 




gh^vv* = -^giy/mh^m-pm-p, , gn^vv = 251 



'Xco 1 _ 



where fy^^g and are the decay constants of Xco and J/V', respectively. ~ 0.41 GeV, which is determined by the 
experimental data [4]. And „ can be approximately determined by the QCD sum rule approach, i.e. ^ 0.51 GeV 

. . \ 

There is no much information on the couplings of a charmed baryon to a charmed meson and light baryon. If 
considering SU(4) flavor symmetry, it will relate different couplings with each other. For the considered = 1/2+ 
charmed baryons, we would expect the following relations [3ll |: 

[t--D+p] = [t--Dt^+] = -V2[S^°Z?+S-] = -V2[%-Dt^'] 
v3 

[A-D+n] = = = -V6[S-i?+A] 

^;[KDtk] = -[pK+k], (19) 

where the square bracket "[• • •]" denotes the coupling constant for the corresponding vertex, and the pseudoscalar 
meson in the bracket can be replaced by the corresponding vector meson. These formulas relate the charmed-baryon- 
meson couplings with the strange-baryon- meson coupling. For instance, the commonly adopted value for AKN 
coupling is g^KN — —13.2 |32h34| |. from which we obtain gA^DN = —givKN — 13.2. Such a relation may contain 
rather large uncertainties as we can see that the QCD sum rule approach suggests a smaller value IgAc-Dwl — 6.7 ±2.1 
[ssj . Experimental data from the charmed meson photoproduction and charmonium absorption by nucleons may 
offer some constraints on these couplings. Unfortunately, the extractions of the couplings still depend on the adopted 
theoretical models [l^, H^H^]. As a tentative solution for this, we empirically retain the SU(4) symmetry relations 
in the numerical calculation. The relevant strange-baryon-meson couplings are taken from the Nigmegen potential 
model as follows [H, : 



9AK'N = -4.26, KAK'N = 2.16, 

9j:k'N = -2.5, K^K'N = -0.22, 



(20) 
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BR(in units of 10"") 


PP 


AA 


S-S+ 




Hadron loop 


9.0 ~ 17.0 


6.3 ~ 12.5 


5.05 ~ 10.0 


4.82 ~ 9.56 


Exp. 


13 ±4 


10.4 ±3.1 







TABLE I: Branching ratios for i]c — )■ YY predicted by the intermediate charmed hadron loop transitions in the range a = 
0.47 ~ 0.53 which corresponds to the measured lower and upper bound of BR(r;c pp) ■ The available experimental data are 
taken from Ref. _4j, and the dashes mean that the data are unavailable. 



and 

9AKN = -13.2, g^KN = 3.9. (21) 

The form factor parameter a generally cannot be determined from the first principle. It will depend on a particular 
process and its value is order of unity. There might exist some differences between the values of the form factor 
parameter for rjc pP and Xco — ^ Pp since r]c and Xco belong to different spin multiplets. It is natural to anticipate 
that the counter term structures would be different for loops involving different spin multiplets. Similar feature was 
also found in Ref. [28| and other studies Q. As emphasized in the literature, in order to reduce the uncertainties 
arising from the form factor, one first has to rely on the experimental data to determine the form factor parameter, 
and then apply it to other flavor-symmetry-related processes to make predictions which can be examined by further 
experimental data. This empirical treatment will also be useful for a further control of the uncertainties arising from 
the SU(4) flavor symmetry breaking. Namely, with the charmed-baryon-meson couplings fixed by the SU(4) relations, 
we allow the experimental data for rjc — ?> pp to determine the range of a as a compensation of the uncertainties from 
the couplings. This can be regarded reasonable since a physical coupling form factor will generally be correlated with 
these two aspects. Reliability of this treatment can be tested by the prediction for the 77c — ^ AA branching ratio, 
which turns out to be consistent with the experimental data 

At present, the branching ratios for r/c — > pp and XcO ~^ Pp have been measured by experiment Q. Signals of he 
were also found in pp annihilations (sgI . [37| , which is a hint that he — ^ pp could be an important channel in the he 
decays. Thus, we will use the measured branching ratios BR(?7c — > pp) and BR(xcO ~^ Pp) to extract the form factor 
parameters for the S'-wave state rjc and P-wave state XcO: respectively. Considering that he belongs to the same spin 
multiplet as XcO; we conjecture that they may share the same intrinsic dynamics in their decays into pp. We will then 
adopt the same form factor parameter for he as that extracted from XcO — >■ pp - 

We use the software package LoopTools to calculate the loop integrals [23 • The results are displayed in Tables HI 
im and mil The experimental data for rje — >■ pp and XcO Pp are adopted for the determination of the form factor 
parameter a with the range of the uncertainties. In Table HI the predicted branching ratio, BR(?7c — > AA) = (6.3 ~ 
12.5) X 10~4, is consistent with the data BR'^^'p (rje -)- AA) = (10.4 ±3.1) x 10"^ which is a sign for the parameter 
under control. The experimental data for rjc — >■ and SS are unavailable. Our calculations suggest that these two 
branching ratios are compatible with that for rje — )■ A A. This expectation can be examined by the BESIII experiment. 

The experimental data for XcO ~^ Pp and AA will allow a further check of the model and its parameter space in 
the spin-1 multiplets. In Table HIl the calculation results are given by the form factor parameter within a range, i.e. 
a = 0.23 ^ 0.24, which corresponds to the experimental uncertainties of XcO ^ Pp [4]. This value range is different 
from that for the rje decays with a = 0.47 ^ 0.53. This can be understood as a consequence that XcO and rjc belong to 
different spin multiplets. We shall show later the sensitivities of the calculation results to the form factor parameter 
later. For the XcO decays, the data from CLEO suggest relatively larger branching ratios for XcO AA compared 
with that for XcO — ^ Pp- Also, branching ratios for XcO — ^ and SS are sizeable. In contrast, with the same form 
factor parameter a, we find relatively smaller branching ratios for XcO — > AA, SE and SS. It could be a sign that 
the SU(4) flavor symmetry is badly broken. Namely, Eq. (IT51) may be too rough, and can only provide an estimate 
of magnitude orders for XcO and S2. In this sense, the calculation results for XcO AA, EE and SS, though 

turn out to be smaller than the experimental data, can be regarded as reasonable. 

Although the present model uncertainties do not allow us to conclude the magnitudes of the hadron loop contri- 
butions, the pattern predicted by this mechanism still suggests an important role played by the hadron loops in the 
explanation of the helicity-selection-rule violations in these exclusive decays. Meanwhile, we stress that it would be 
essential to have improved experimental measurements in order to gain better insights into the transition mechanism. 
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BR(in units of 10"'*) 


PP 


AA 


E-E+ 




Hadron loop 


1.96 ~ 2.34 


1.19 ~ 1.51 


0.55 ~ 0.69 


0.52 ~ 0.66 


Exp. \£ 


2.15 ±0.19 


4.4 ± 1.5 




< 10.3 


Exp. [39| 


2.25 ±0.27 


4.7 ± 1.6 


3.25 ± 1.14 


5.14 ± 1.25 



TABLE II: Branching ratios for XcO YY predicted by the intermediate charmed hadron loop transitions in the range 
a — 0.23 ~ 0.24 which corresponds to the measured lower and upper bound of BR(xcO ^ PP) 4]. 



BR(in units of 10"") 


PP 


AA 


E-S+ 




Hadron loop 


15.2 ~ 19.3 


5.88 ~ 7.47 


4.56 ~ 5.80 


5.57 ~ 7.08 


Exp. 











TABLE III: Branching ratios for he — >■ YY predicted by the intermediate charmed hadron loop transitions. The a range for 
he — >■ YY is taken the same as that for XcO YY, i.e. a = 0.23 ~ 0.24. The dashes mean that the data are unavailable. We 
take the width of he as T{hc) = 0.73 MeV, which is the central value measured by BESIII recently [iOj . 

Experimental data for he — >■ YY so far are unavailable. Thus, our predictions are based on the assumption that 
he shares the same form factor parameter as XcO since they belong to the same spin multiplet. In Table IIIII we 
list the branching ratios given by a = 0.23 ~ 0-24, which is the same as adopted in XcO PP- In comparison 
with the theoretical calculations of Refs. [l^, Hi], [42], our prediction of the branching ratio of he — )■ pp seems to 
be larger. In fact, the theoretical predictions in the literature also appear to be quite different from each other. 
In particular, some of those results strongly depend on the evaluation of BR{hc J/V'tt") in the combined cross 
sections for pp ^ he ^ J/ i/)7r° from the E760 data [H, [H, [4l| . Another reason for the discrepancies among these 
theoretical estimates may be due to different intrinsic mechanisms adopted for the explanation of the helicity-selection- 
rule violation. We expect that the future precise measurement of he — > YY will help disentangle the underlying 
mechanisms. 

In Figs. [51 E]and[7J we also examine the dependence of the results on the form factor parameter a. The adopted 
values of a are within a reasonable range and well controlled by the available experimental data. Although some 
uncertainties will be inevitably introduced by the phenomenological form factor, the branching ratio fractions among 
the considered channels turn out to be stable and less model-dependent. This feature suggests that the branching 
ratio fractions are less model-dependent in comparison with the absolute branching ratios. In another word, although 
the model predictions for the absolute branching ratios are lack of experimental constraints, thus, becomes sensitive 
to the form factor parameter a, we would expect that the predicted branching ratio fractions among those considered 
channels are less sensitive to it. As a consequence, any experimental results for BR(/ic YY) at 0(10"'') ^ 0(10"'^) 
would imply the importance of the hadron loop contributions to this helicity-selection-rule violation transition. 

In the P-wave charmonium decays the next higher Fock state ccg, i.e. the so-called color octet, will contribute 
at the same order as the color singlet cc in the framework of perturbative factorization method |43| - |46| . This 
scenario may share the same intrinsic physics with those intermediate charmed hadron loop transitions based on 
quark- hadron duality argument, i.e. a manifestation of the same physics at either quark- giuon level, or hadron level. 
Some qualitative discussions can be found in Refs. [13, H^, 113) and references therein. We address that since it is not 
easy to handle these exclusive decays with the perturbative methods considering the relatively low energy scale, the 
intermediate charmed hadron loop transition would serve as a natural soft mechanism in the numerical exploration. 

IV. SUMMARY 

In this work, with an effective Lagrangian method based on heavy quark and flavor symmetry, we investigate the 
role played by the intermediate charmed hadron loops in the processes rje, XcO, he — 5> YY , which are supposed to be 
highly suppressed by the helicity selection rule. The results indicate that the transitions via these kinds of loops 
as long-distance effects can give significant contributions. This is a further test of the mechanism for the evasion of 
helicity selection rule that we proposed in Ref. iQJ . Although the model bares uncertainties arising from the unknown 
coupling constants and form factor parameter, the available data have provided a reasonable control on the range of 
the form factor parameter values. Branching ratios for some unmeasured channels can thus be predicted. Sizeable data 
samples on charmonium decays accumulated at BEPCII and CLEO-c, and future proton-antiproton annihilation data 
from Panda, are expected to provide a great opportunity for revealing the underlying mechanisms for charmonium 
hclicity-selection-rule-evading decays. 
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a 



FIG. 5: a-dependence of the calculated branching ratios for rjc — >■ YY. 




a 



FIG. 6: a-dependence of the calculated branching ratios for XcO — > YY . 
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